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ABSTRACT 



We report the results of a numerical search for non-topological 
solitons in the two-Higgs standard model, characterized by the 
non-trivial winding, 7r3(S'^), of the relative phase of the two dou- 
blets. In a region of (weak-coupling) parameter space we identify 
a branch of winding solutions, which are lower in energy than 
the embedded standard sphaleron and deformed (bi)sphalerons. 
Contrary, however, to what happens in 2d toy models, these so- 
lutions remain unstable even for very large Higgs masses. 



Two of us have analyzed recently non-topological winding solitons in 
renormalizable gauge models with two Higgs scalars. Such solitons arise in 
two and three space-time dimensions [0], and can be embedded as extended 
membrane defects in the weakly-coupled two-Higgs standard model (2HSM) 
[0]. It was furthermore argued by analogy ^ that the 2HSM may have string- 
like or particle-like excitations characterized by a non-trivial mapping of the 
spatial three sphere onto the relative phase of the two doublets. Depending 
on the details of the potential, the relevant quantity could be either the Hopf 
invariant, 713(6'^), or the winding number, n^lS^). Here we report the results 
of a numerical search for solitons of the latter type. The search revealed 
indeed a branch of winding solutions, but unlike their 2d analogs these stay 
classically unstable for arbitrarily large Higgs masses. These new solutions 
could nevertheless be of interest for the study of electroweak baryogenesis in 
the 2HSM p. Sphaleron solutions in the 2HSM have been discussed previ- 
ously in ref. ^ and in a different context in ref. 0. 

There exists a close analogy between sphaleron solutions in the minimal 
standard model 0; HI; and those in the abelian-Higgs [|1^] or global double- 
well models on the circle. Similarly, the structure of classical solutions in 
the 2HSM is strongly reminiscent of that in the abelian two-Higgs, or in the 
global mexican-hat model on the circle. This latter model is defined by the 
action 

, (1) 

i2 ^ 4^ ^ J ' ^ ^ 

whose single classically-relevant parameter is mL = V2XvL. We henceforth 
set L = 1. Static solutions correspond to motions of a particle in the in- 
verted mexican-hat potential with period 27r. One solution, which exists for 
all values of m and can be thought of as the basic sphaleron [|, corresponds to 
the particle sitting at the bottom of the inverted hat. It has energy 7rm^/8A, 
starts out with two negative modes, and acquires two extra ones every time 
m/\/2 crosses a positive integer n. Two extra branches of solutions bifurcate 



r2TTL 

S = I dt dx 



* Though it does not in this case represent a minimum-height passage between two 
discrete degenerate vacua. 
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at these points: the Sn branch describes an n-fold oscillation about the bot- 
tom of the inverted hat, with zero angular momentum. These solutions, also 
present in the double-well model, are everywhere unstable and are the analogs 
of the deformed sphalerons of the standard model 0. The Wn branches 
on the other hand describe an n-fold rotation around the tip of the hat, and 
have no analog in the one-Higgs model. The solutions have a constant scalar 
magnitude = (m^ — 2n^)/2\, an energy equal to 7m'^{m'^ — n^)/2\, and 
become classically stable beyond the critical values = Qri^ — 1 New 
branches, Wn-, emerge at these points; they have slightly higher energy and 
correspond to winding motions with an oscillating scalar magnitude. Further 
bifurcations occur along the Wn>i branches, and possibly also along Sn and 
Wn, but an exhaustive analysis is beyond the scope of this letter. Figure 1 is 
an illustration of the topological tree of solutions. The main lesson to retain 
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Figure 1. Part of the tree of classical solutions for the model (1). The vertical axis gives 
the energy, measured in units of m-^/SA so as to simplify comparison with ref. The W\_ 

branch is drawn out of scale: its energy differs by only ~1% from the energy of Wi. 
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from this toy model is that for m > \pl the winding solution on the W\ 
branch is the one with lowest energy, and that furthermore for m > -\/5 it 
becomes a classically-stable soliton. 

Using this toy model as a guide, we turn now to the 2HSM. We will 
work with the simplifying assumptions of vanishing Weinberg angle and of a 
potential consisting of two independent mexican hats. We believe that these 
restrictions are optimal, in the sense of favouring the stability of the sought- 
for winding solutions. We also restricted our search to spherically-symmetric 
configurations. The action is 



S 



^ a=l,2 



A 



2\2 



a=1.2 
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(2) 



1,2, and Wf,u 



with D^Ha = {d^ + W^)Ha, for a 
the length scale by setting mw = g\l {vf + v^)!'! 
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We fix 
and scale out of 



the action the semiclassical parameter aw = 5'^/4vr. There remain three 
classically-relevant parameters: the two radial Higgs masses rua = 2i/A^i'aj 
and tan/5 = V2/V1. 

The general spherically-symmetric static ansatz, following the conven- 
tions of Yaffe IHl, reads 



Wi = - 



1 

2i 



Wq = woTini/2i , 
(a - l)eijkrjnk/r + h{6ij - ninj)Tj/r + WiUirijTj 

Ha = (/ia + il^ariiT-i)^ 



(3) 



where nj is the unit vector in the direction of r, Xj are the three Pauli 
matrices, ^ is a constant unit doublet and wq, wi, a.,h, fia,J^a are functions 
of r to be determined by the stationarity of the energy. It is convenient to 
choose the gauge Wq = 0, and then solve Gauss' constraint to obtain wi = 
p. Defining the complex fields (pa = fJ^a + i^a = /aGxp {i9a) (a = 1,2) and 
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X = a + ib = Ixl exp (iip), one can put the energy functional in the form: 

1 /"oo f 1 

E =— / dr\ \xf + 2rW/3(|0;p + \<P',\') + ^,{\x\' - if 
aw Jo [ 

+ cos2/5(|xr + l)(|0ir + 1021') - 2cos2/? /2e[x*(0? + 0^)] (4) 
+ \mycos^mi? - I? + ^^2r'^(l02p - tan2/3)2|, 

with primes standing for derivatives with respect to r. The above expression 
is invariant under a simultaneous rotation of 0i, 02 and x by a phase uo^uj 
and 2u; respectively, as well as under independent changes of sign of 0i and 
02- These symmetries together with the requirement of finite energy allows 
us to fix the boundary conditions at r = cxd in the form: 

X — >■ 1 , 01 — 1 and 02 — > tan/3 . (5) 

At r = on the other hand, finiteness of the energy, of the gauge current 
and of the field strength imply |]^: 

Ixl ^ 1 , X*0a ^ I0a|' and x' ^ . (6) 

Notice in particular that at the origin either fa = 0, or else 26a — = 
(mod 27r). For configurations with non- vanishing Higgs magnitudes we 
may conclude that 

61-62 = Ntt (7) 

r=0 

The integer N is the gauge-invariant winding number characterizing the map- 
ping from the spatial onto the SU(2) manifold of the relative phase of the 
two doublets 

We used the relaxation method [|rT|] to integrate the field equations nu- 
merically. Starting with some initial guess {$} for the solution, one linearizes 
the static field equations, solves them, and adds the result to {$} thus ob- 
taining an improved guess for the next iteration. The radial coordinate r 
was replaced by n-(-l=201 points r^, with tq = and Vn = R = 20. This was 
sufficient to achieve an accuracy of order 10^^ to 10~^ in the energy of the 
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solutions. Following ref. |Q we distributed the lattice points according to 
the formula = (1 - Q^s/ny (^^ks/n _ ^-^^^ ^ ^ y^^f^^ ^ fiR)/{l + R)] and 
fi = max(mi,m2). The linearized equations at each iteration step take the 
form Efj^*^ = e}^^ where 5$ is a vector of dimension n x (^fields) = 1200, 
normalized so that the kinetic energy reads ^Y.^^^ and E^^\ E^'^^ are 
respectively the first and second variations of the static energy, eq. (4). The 
matrix E^"^^ has a special block-three-diagonal form, each block correspond- 
ing to the six fields at the same lattice point. We have exploited this special 
form to invert the matrix by a forward-elimination and back-substitution 
procedure. This same procedure also allowed us to calculate the spectrum of 
negative fluctuation modes around the solutions. 

Figure 2 summarizes schematically the tree of solutions that we have 
found in the particular case of equal Higgs masses, mi = m2 = m. This 
further restriction of parameter space has the following convenient feature: 
any solution {H, Wfj_} , of the one-doublet model with Higgs mass m and 
niw = 1, yields a solution {Hi = Hcos(3, H2 = Hsinf], = W^} of the 
2HSM with the same total energy though a different fluctuation spectrum. 
The branches Sq and 5*1 in particular correspond to the standard sphaleron 
[^] and deformed (bi)sphalerons ^ of the minimal model, and have been 
used to check the accuracy of our numerical routines. All other branches 
shown in figure 2 have no analog in the one-doublet model. 

The standard sphaleron (with x, (pa real) is the only solution we have 
found for m/my/ < 5.5. It has a single mode of instability in this range, but 
develops a second one where the new branches of winding solutions, {Wi) and 
its conjugate, bifurcate. Extra negative modes appear at each subsequent 
bifurcation: at m/mw — 12 (first deformed sphalerons), at m/mw — 50 
(doubly- winding solutions), at m/mw — 138 (second deformed sphalerons) 
and so on down the line. This sequence of bifurcations is strongly reminiscent 
of the 2d toy model, except that the Wn and Sn branches do not emerge in this 
case simultaneously. Notice that the topological invariant, solutions 
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Figure 2. Solutions of the two-Higgs model, eq. (2), for tanl3=l, as function of the 
(common) radial Higgs mass m. Both the vertical and the horizontal axes have non-uniform 
scales. The number of negative modes is shown explicitly for every branch. The winding branch, 
Wi, emerges out of the main sphaleron at m/niw—^-^-, but remains unstable at least up to a 
mass ratio of ^^200. 

with Us tlie number of unstable modes of tlie solution, is conserved at all 
bifurcation points. This is also true at m/mw — 14.9, where a once- and 
a twice-unstable branch emerge out of the first deformed sphaleron, corre- 
sponding to the transformation of a local maximum, into two maxima sepa- 
rated by a lower-energy saddle. 

The most interesting feature of Figure 2 is of course the winding branch 
Wi, which on the basis of the 2d analog could have evolved into a stable 
soliton at sufficiently large Higgs masses. The negative eigenvalue does not, 
however, cross the zero axis in this case, at least all the way out to m/mw ~ 
200 (see table 1) f\. By performing a numerical scan, we have checked that 

^ Such large values of parameters are in any case academic, since the semiclassical 
approximation cannot be anymore trusted. 
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the winding branch does not stabihze anywhere inside the larger space of 
parameters of the potential, eq. (2) (see Table 2). Since we do not expect 
interactions between the two doublets to improve stability, we beheve that 
such winding solitons do not exist. This is the main conclusion of the present 
letter. 

The unstable solutions are nevertheless of interest in their own right. 
They emerge inside the region of weak coupling, and can be lower in energy 
and comparatively less unstable than the known sphaleron and deformed 
sphaleron solutions, as shown in Tabic 1. They furthermore continue to exist 
in a large region of parameter space, which is sampled selectively in table 2. 
The profile of a typical solution for tan/5 = 1.4, mi = 0.1 and 1712 = 10.1 
is plotted in Fig. 3. As is the case with deformed sphalerons in the mini- 
mal model, the appearance of the two conjugate Wi branches corresponds to 
the splitting of a minimum-height passage in the energy landscape into two. 
These configurations should thus dominate baryon-number violating transi- 
tions at high temperature, in the appropriate regions of parameter space of 
the 2HSM. More generally, their existence illustrates how rich the configu- 
ration space of the 2HSM is. Its exploration could, we believe, reveal many 
surprising features. 



mi = m2 
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6.09 


30.0 
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4.75 


2.42 


5.25 


291.8 


5.06 


4.97 


100.0 


4.75 


2.41 


5.33 


1192 


5.07 


4.75 



Table 1. The energy (in units of mw/ciw) and negative eigenvalue of the fluctuation 
spectrum (in units of mw) along the winding branch Wi, in the symmetric case tanf3=l and 
mi=m2. We have included for comparison the energy and most negative eigenvalue along the 
embedded sphaleron and first deformed sphaleron branches, respectively So and Si. 
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mi 


m2 


tanj3 






10.0 


10.0 


0.03 


4.78 


11.20 


10.0 


10.0 


0.1 


4.78 


11.14 


10.0 


10.0 


0.5 


4.74 


5.43 


10.0 


10.0 


1.0 


4.67 


3.09 


10.0 


10.0 


1.4 


4.70 


3.70 


10.0 


10.0 


1.8 


4.73 


4.82 


10.0 


10.0 
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4.53 
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4.51 
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4.49 


5.94 


3.0 


6.0 
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4.46 


5.85 


2.5 


6.0 


2.0 


4.44 


5.80 


2.0 


6.0 


2.0 


4.43 


5.76 


1.5 


6.0 


2.0 


4.40 


5.72 


10.0 


0.0 


4.0 


3.17 


1.42 


10.0 


0.1 


4.0 


3.26 


1.58 


10.0 


0.4 


4.0 


3.45 


1.92 


10.0 


1.0 


4.0 


3.71 


2.45 


10.0 


2.0 


4.0 


4.00 


3.21 


10.0 


5.0 


4.0 


4.45 


5.29 


10.0 


10.0 


4.0 


4.77 


9.83 


10.0 


15.0 


4.0 


4.90 


6.43 


10.0 


20.0 


4.0 


4.95 


5.25 


10.0 


30.0 


4.0 


4.99 


4.67 


10.0 


10.0 


1.0 


1.99 


1.51 



Table 2. The energy and negative eigenvalue of the W\ branch, for various values of 
tan/3 and of the ratio m2/mi. The units are the same as in table 1. Notice that when tanP—>0 
the winding solution merges with the sphaleron. For equal values of the two Higgs masses the 
lowest-energy, least unstable solution is obtained in the symmetric case tan/3=l, when the winding 
is equally shared by the two Higgses. 
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Figure 3. The profile of a solution on the Wi branch for Higgs masses about 8 and 800 
GeV respectively. 
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